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1.  Introduction 

Safety  margins  protect  engineering  devices  against  extreme 
load  events,  unexpected  operating  conditions,  and  manufacturing 
imperfections.  Early  design  rarely  takes  into  account  reliability 
and,  as  a  result,  safety  factors  are  applied  as  a  posteriori  correc¬ 
tions;  these  are  typically  based  on  experience  and  change  radi¬ 
cally  In  different  application  areas. 

There  is  a  general  lack  of  well-defined  procedures  to  objec¬ 
tively  quantify  the  confidence  in  the  operability  of  an  engineering 
system.  Complex  devices,  such  as  aerospace  systems,  are  typically 
difficult  to  characterize  in  a  satisfactory  way.  The  behavior  of  a 
complex  engineering  device  is  also  typically  characterized  by  a 
large  number  of  parameters;  failures  typically  occur  in  response 
to  one  of  the  operating  conditions  exceeding  a  permissible 
threshold,  e.g.  exceeding  structural  load.  In  addition,  unexpected 
(hidden)  correlations  between  physical  processes  can  lead  to 
failure  in  response  to  a  sequence  of  events  with  each  parameter 
remaining  within  the  allowable  range  of  operation. 

The  operation  region  is  the  environment  in  which  the  device  of 
interest  is  designed  and  expected  to  operate  successfully;  math¬ 
ematically  we  define  this  region  as  the  admissible  volume  Qor  in 
the  space  of  all  the  parameters  affecting  the  behavior  of  the 
system.  For  is  the  boundary  of  such  a  region  and  is  defined  as  the 
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locus  of  the  unsafe  operating  conditions;  these  do  not  necessarily 
represent  failure  points,  but  cases  for  which  no  positive  evidence 
of  acceptable  performance  exists;  we  refer  to  them  more  precisely 
as  the  performance  thresholds. 

The  Intuitive  notion  of  safe  operation  of  the  device  implies  a 
distance  (margin)  from  the  performance  threshold  and  therefore  a 
parameter  space  Q;  smaller  than  the  operating  region  (£2s  c  Qor, 
see  Fig.  1 ).  The  main  difficulty  In  determining  the  margin  Is  the 
presence  of  uncertainties  in  both  the  precise  location  of  the  cliffs 
(For)  and  the  actual  operating  conditions  of  the  system  (the 
location  in  Qor).  These  uncertainties  are  due  to  variability  in  the 
environment  the  system  is  operating  in,  manufacturing  toler¬ 
ances,  material  imperfections  {aleatory  uncertainties),  and  to  our 
limited  understanding  of  the  physical  processes  at  play  (epistemic 
uncertainties).  As  a  result,  a  purely  deterministic  evaluation  of  the 
margin  is  not  sufficient  and  probabilistic  approaches,  or  more 
generally,  approaches  that  directly  account  for  the  uncertainties 
have  to  be  considered. 

Quantification  of  Margin  and  Uncertainties  (QMU)  is  a  meth¬ 
odology  created  to  facilitate  analysis  and  communication  of 
confidence  for  certification  of  complex  systems  [1].  In  QMU  the 
confidence  is  defined  in  a  deceptively  simple  way: 


where  M  Is  the  a  measure  of  the  margin  and  U  a  measure  of  the 
uncertainty.  CR  is  the  confidence  ratio  which  has  to  be  evaluated 
for  the  full  system  (and  for  the  subsystems)  in  the  operating 
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region  Qs  of  the  device.  A  CR  sufficiently  larger  than  one 
intuitively  indicates  safe  conditions. 

The  theoretical  framework  for  QMU  and  the  open  questions 
related  to  Its  applicability  are  discussed  in  a  number  of  articles 
and  reports,  for  example  [2-6].  Detailed  applications  of  QMU  to 
real-world  engineering  studies  are  limited  because  of  the  diffi¬ 
culties  in  assessing  the  predictive  capabilities  of  complex  compu¬ 
tational  models.  In  [2]  several  applications  ranging  from  diffusion 
in  porous  media  to  interactions  and  reflections  of  shock  waves  in 
detonation  problems  are  presented  together  with  the  analysis  of 
the  importance  and  the  potential  impact  of  errors. 

A  rigorous  methodology  to  study  failure  analysis  is  introduced 
in  [6]  for  a  simple  problem  of  imploding  and  exploring  ring 
structure.  The  objective  is  to  demonstrate  that  the  system 
performance  under  a  worst-case  scenario  remains  above  a  given 
threshold.  This  QMU  framework  is  used  to  explore  different 
situations  in  which  the  computational  model  Is  assumed  to  be 
either  exact  or  imperfect  and  various  scenarios  for  uncertainty  in 
the  performance  measures  are  considered.  In  particular,  the 
authors  explore  how  experimental  data  can  be  used  to  comple¬ 
ment  non-exact  simulations.  The  QMU  framework  is  based  on 
formal  bounds  for  the  probability  of  exceeding  stresses  on  the 
structure. 

A  more  complex  QMU  analysis  is  reported  in  [5]  where  the 
thermal  response  of  an  engineering  device  under  fire  is  studied. 
The  complexity  of  the  simulations  and  the  need  to  rely  on  low 
fidelity  computational  models  representing  the  system  Is  dis¬ 
cussed  together  with  the  importance  of  calibration  and  the 


Fig.  1.  Regions  depicting  the  parameter  space  of  an  engineering  system.  £3„r  is  the 
operating  region,  For  represents  the  performance  threshold  and  Qs  the  safe 
operation  region. 


difficulty  in  computing  tails  of  probability  distributions  (failure 
probability). 

The  objective  of  this  paper  is  to  illustrate  how  a  QMU  analysis 
is  formulated  in  a  complex  engineering  problem  where  the 
interplay  between  fluid  dynamics  and  thermal  loads  determine 
a  performance  threshold.  In  the  following  we  will  first  briefly 
describe  the  operability  limits  of  an  air-breathing  hypersonic 
vehicle  and  then  Introduce  a  relevant  QMU  framework.  We  will 
then  describe  a  simple  one-dimensional  computational  model  of 
the  propulsion  system  In  order  to  explore  the  uncertainties 
associated  with  the  system  predictions.  Finally,  we  will  determine 
the  safe  operating  region  for  a  range  of  fuel  flow  rates  and 
combustor  configurations. 


2.  Operability  limits  of  an  air-breathing  hypersonic  vehicle 

Air-breathing,  hypersonic  vehicles  are  highly  Integrated  sys¬ 
tems  whose  performance  depends,  to  a  very  large  extent,  on  the 
complex  physics  and  interactions  between  all  of  their  compo¬ 
nents.  Such  performance-critical  systems  cannot  be  credibly 
designed  with  today’s  state-of-the-art  simulation  capabilities 
not  only  because  there  are  no  available  methods  to  predict  the 
basic  physical  phenomena,  but  also  because  of  uncertainties 
present  in  the  actual  flight  environment,  in  the  manufacturing 
of  the  various  components,  etc.  As  a  result,  the  final  design  relies 
on  extensive  physical  testing  [7]  and  the  Introduction  of  con¬ 
servative  assumptions  and  safety  constraints. 

Of  particular  interest  in  this  paper  are  the  operability  limits 
associated  with  the  propulsion  system  of  hypersonic  vehicles 
flying  at  mid-range  Mach  numbers  (Mass  8).  In  these  circum¬ 
stances,  the  engine  operates  under  supersonic  conditions  - 
scramjet  mode  -  and  a  key  performance  metric  is  the  amount  of 
heat  released  In  the  combustion  chamber,  as  it  is  directly 
connected  to  the  generation  of  thrust. 

Fig.  2  shows  a  schematic  of  the  cross-section  of  a  hypersonic 
vehicle.  The  presence  of  a  supersonic  flow  stream  over  the  entire 
vehicle  logically  Introduces  a  streamwise  splitting  of  the  vehicle 
into  several  subsystems. 

The  forebody  Is  dominated  by  the  presence  of  the  strong  bow 
shock  and  the  associated  complex  thermo-chemistry  effects.  The 
extreme  heating  occurring  at  the  stagnation  point  determines  the 
choice  of  materials  and  cooling  strategies.  As  the  flow  decelerates 
towards  the  engine,  boundary  layers  develop  on  the  vehicle 
surface.  Turbulence  trips  are  designed  to  force  transition  to 
turbulence:  in  spite  of  the  increased  surface  heating  and  friction, 
turbulent  boundary  layers  entering  the  engine  Inlet  lead  to 
increased  stability  and  mixing.  The  inlet/isolator  system  is 
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Fig.  2.  Schematic  representation  of  an  air-breathing  hypersonic  vehicle. 
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designed  to  condition  the  flow  towards  pressure  and  temperature 
conditions  that  are  most  favorable  to  combustion.  In  addition, 
the  shock  train  present  in  the  isolator  reduces  the  distortion  of 
the  incoming  flow  due  to  angle  of  attack  and  yaw  variability.  The 
combustor  is  the  most  critical  component  of  the  vehicle  and,  in 
particular,  the  injection  system  is  carefully  designed  to  mix  the 
fuel  and  the  incoming  supersonic  air  stream  and  to  produce  auto¬ 
ignition  of  the  mixture.  Finally,  the  nozzle  and  afterbody  further 
accelerate  the  flow  stream  producing  thrust. 

One  of  the  main  challenges  in  designing  a  scramjet  engine  is 
the  extremely  short  residence  time  of  the  fuel  within  the 
combustor  ( —  ms);  the  ability  to  inject  the  fuel,  mix  it  in  the 
supersonic  incoming  air  stream  and  ignite  the  resulting  mixture  is 
the  critical  design  objective  for  the  propulsion  system.  A  mini¬ 
mum  combustion  efficiency  (or  heat  release)  is  a  performance 
threshold  for  the  vehicle  design.  However,  there  is  a  well- 
identified  limit  to  the  heat  release  that  can  be  deposited  in  a 
supersonic  stream  [8].  Excessive  heating  leads  to  thermal  choking: 
a  normal  shock  and  a  consequent  subsonic  flow  region  is 
established  in  the  combustor.  In  addition  to  a  reduction  in  the 
performance,  this  can  lead  to  increased  structural  and  thermal 
loads  and  eventually  lead  to  failure.  The  normal  shock  can  also 
propagate  upstream  in  the  combustion  chamber,  and  eventually 
interact  with  the  isolator,  creating  extensive  regions  of  boundary 
layer  separation;  moreover,  the  shock  motion  can  lead  to  engine 
unstart  with  the  entire  isolator  shock  train  moving  upstream. 
Under  these  conditions  the  vehicle  performance  is  compromised 
and  extreme  actions  have  to  be  taken  to  restart  the  engine  [7]. 

Unstart  conditions  can  also  be  reached  for  different  reasons, 
not  directly  connected  to  the  combustion  process.  In  particular, 
indications  of  unstart  events  connected  to  perturbations  of  the 
flow  at  the  inlet  [9]  and  to  thermal  deformations  of  the  structure 
[10]  have  been  observed  in  ground  tests. 

The  objective  of  the  present  analysis  is  to  characterize  the 
thermal  choking  conditions  and  to  quantify  the  margin  associated 
with  the  unstart  limit.  In  particular,  it  is  clear  from  the  above 
description  that  the  engine  performance  (measured  in  terms  of 
thrust)  increases  with  the  heat  addition  until  it  reaches  the 
critical  unstart  limit,  so  it  is  desirable  to  operate  the  engine  as 
close  as  possible  to  this  cliff. 

In  the  following  we  will  introduce  a  highly  simplified  compu¬ 
tational  model  to  study  scramjet  performance  and  then  we  will 
describe  the  QMU  framework  for  studying  the  confidence  in  the 
predictions. 


2.3.  Sources  of  uncertainty  in  the  simulation  of  vehicle  unstart 

The  unstart  phenomenon  is  a  complex  thermal/fluid  dynamic 
process  that  is  critically  sensitive  to  a  variety  of  geometrical  and 
physical  parameters.  It  is  useful  to  provide  a  characterization  of 
the  major  sources  of  uncertainties  in  this  process  to  guide  the 
development  of  the  QMU  framework. 

We  can  distinguish  between  various  sources  of  uncertainty.  The 
first  group  includes  the  imprecise  characterization  of  the  environ¬ 
ment  in  which  the  vehicle  is  flying.  Of  these,  we  can  assume  that  to 
leading  order,  the  speed  of  flight,  the  angle  of  attack  and  the  overall 
atmospheric  conditions  (mainly  in  terms  of  temperature  fluctua¬ 
tions)  will  have  the  highest  impact  on  the  unstart  predictions. 
Another  large  set  of  uncertain  parameters  is  related  to  the  material 
and  static  characterization  of  the  vehicle.  In  this  category  we  can 
include  structural  and  material  inhomogeneity,  fuel  mixture 
imperfections,  surface  roughness  and,  more  generally,  out-of-spec 
geometry.  The  final  set  of  parameters  are  related  to  the  thermal 
and  fluid  processes  during  flight:  fuel  injection  rate,  fuel  tempera¬ 
ture,  combustion  process,  thermodynamic  non-equilibrium  at  the 


leading  edge,  and  the  laminar-to-turbulent  transition  at  the 
engine  inlet. 

Numerical  predictions  of  hypersonic  vehicle  operations  are 
based  on  computational  tools  that  require  models  to  simulate 
physical  processes  that  cannot  be  treated  from  first  principles. 
Therefore,  a  considerable  level  of  uncertainty  is  also  associated 
with  the  hypotheses  and  assumptions  used  in  such  models. 

3.  A  reduced-order  model  for  computing  scramjet 
operability  limits 

Engineering  analyses  of  scramjet  propulsion  systems  have 
been  carried  out  in  the  past  [11]  with  the  objective  of  designing 
access-to-space  vehicles  and  high-speed  commercial  cruise  air¬ 
planes.  A  comprehensive  review  of  design  methodologies  is 
reported  in  Heiser  et  al.  [12]:  semi-analytical  relations  based  on 
compressible  gas-dynamics  equations  are  derived  and  used  to 
illustrate  the  various  design  choices. 

The  objective  of  the  present  analysis  is  not  to  design  a  new 
system  but  to  analyze  how  computational  tools  in  conjunction 
with  experimental  evidence  can  be  used  to  identify  safe  operating 
conditions.  The  need  to  explore  the  thermal  choking  limit 
requires  the  construction  of  a  phenomenological  model  that 
represents  the  unsteady  nature  of  the  flow  in  the  scramjet  and, 
eventually,  unstart  events.  We  seek  the  simplest  approach  that 
can  represent  an  internal  compressible  flow  with  cross-sectional 
area  change,  heat  release,  unsteady  effects  and  the  occurrence  of 
shock  waves.  The  model  we  adopted  is  similar  to  that  introduced 
by  Bussing  and  Murman  [13]  and  represents  only  the  isolator,  the 
combustion  chamber  and  the  nozzle  (see  Fig.  2). 

The  governing  equations  are 

Uf  +  Fx  =  (Z{P— F)  +  Q/Lc  (1) 

with 
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where 


Qcomb  =fstHfmfg(x/Lc)  =  ^fstHfthogix/Lc)  (3) 

p,  u,  e,  p  are  the  local  density,  velocity,  total  energy  and  pressure, 
respectively.  Fx  is  the  convective  flux,  whereas  a(P-F)  on  the  right 
hand  side  represents  the  effect  of  the  cross-sectional  area  change, 
A(x),  with  (Z  =  (l/A)(dA/dx);  Q  is  the  length  of  the  combustor. 
Finally  Q,  represents  the  heat  addition  with  rhf  and  rho  being  the 
mass  flow  rate  of  the  injected  fuel  and  the  incoming  air, 
respectively. /st  is  the  stoichiometric  fuel/air  ratio  and  Ff/the  fuel 
heating  value.  In  the  following  analysis  we  will  consider  hydrogen 
as  the  fuel  and,  therefore,  /st=0.029  and  Hf  =  1.2  x  lO^J/kg,^^). 
The  equivalence  ratio  </>  controls  the  amount  of  fuel  injected. 

The  function  g(xlLc)  in  Eq.  (3)  represents  the  distribution  of 
heat  release  within  the  combustion  chamber  which  is  controlled 
by  the  precise  fuel  injector  geometry,  the  turbulent  mixing,  and 
the  kinetics  of  the  auto-ignition  and  combustion  process.  In  the 
present  one-dimensional  model,  no  attempt  can  be  made  to 
represent  any  of  these  phenomena  realistically  and,  therefore, 
we  assume  a  simple  form  of  the  heat  release  similar  to  [8[: 

g(x/Lc)=[C,^J  (4) 

where  Q  and  C2  are  constants  to  be  determined  from 
experimental  data. 

The  model  is  obviously  a  very  crude  approximation  of  the 
complex  fluid  flow  environment  within  a  scramjet.  It  is  worth 
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mentioning  the  major  assumptions  that  have  been  made  (in 
addition  to  the  one-dimensionality  constraint):  (i)  absence  of 
viscous  effects,  (ii)  adiabatic  walls,  (hi)  no  mass  addition  due  to 
the  fuel  injection  and  (iv)  simplified  heat  release  distribution  due 
to  the  combustion  process.  It  is  possible  to  eliminate  some  of 
these  assumptions;  for  example,  the  original  formulation  in  [13] 
includes  the  effects  of  the  viscous  boundary  layers  on  the  walls  by 
considering  an  effective  cross-section  area;  it  is  beyond  the  scope 
of  the  present  analysis  to  formulate  a  more  comprehensive 
reduced-order  model. 

The  equations  are  solved  using  a  second-order  finite  volume  .M, 
discretization  and  an  implicit  time-integration  scheme.  Convec-  p 

tive  fluxes  are  evaluated  using  the  approximate  Riemann  solver  of  p 

Roe  [14].  Boundary  conditions  are  specified  as  fully  supersonic  at 
the  isolator  Inlet  with  specified  Mach  number,  Ma„  total  pressure, 

Pti,  and  total  temperature  Tt,.  The  outflow  conditions  at  the  exit 
(nozzle)  are  specified  as  either  supersonic  or  subsonic  if  thermal 
choking  occurs. 

In  the  following,  we  describe  the  formal  process  we  followed 
to  verify  and  validate  the  present  reduced-order  model. 

3.1.  Scramjet  geometry  and  design  conditions 

Propulsion  systems  for  hypersonic  air-breathing  vehicles  are 
extremely  simple  systems  compared  to  conventional  jet  engines; 
they  have  no  moving  parts,  simple  geometry  and  a  very  limited 
number  of  parts.  In  particular,  the  configuration  is  typically  a 
straight  or  diverging  duct  (or  pipe)  with  an  end  nozzle  that 
produces  most  of  the  thrust.  The  fuel  injection  system  is  rather 
elaborate  to  improve  mixing;  the  combustor  is  also  a  divergent 
duct  to  enhance  stability  and  prevent  unstart  (we  will  exploit  this 
in  the  following  analysis  of  the  operating  margins). 

The  computational  model  described  above  creates  severe  lim¬ 
itations  in  our  ability  to  represent  the  complex  physical  mechan¬ 
isms  in  a  scramjet  propulsion  system  but  it  allows  us  to  capture 
the  overall  geometrical  configuration;  specifically,  we  consider  the 
geometry  reported  in  [8];  a  constant  cross-section  isolator  of 
length  0.5  m  is  followed  by  a  short  combustor  (Lc  =  0.1  m)  and  a 
nozzle.  Both  the  combustor  and  the  nozzle  are  simple  diverging 
ducts  with  angles  of  Xcomb  =  7.5"  and  a„ozzie  =  f  5  ,  respectively. 

We  will  further  assume  that  the  engine  inflow  conditions 
correspond  to  Ma,=2  at  the  isolator  entrance. 

3.2.  Unstart  scenarios 

In  spite  of  the  simplicity  of  the  present  model,  we  can 
represent  realistic  thermal  choking  scenarios;  in  particular  we 
can  investigate  a  staged  increase  of  the  fuel  flow  rate  within  the 
engine.  The  fuel  addition  leads  to  an  increased  thrust  and  the 
vehicle  operates  within  the  designed  conditions  until  thermal 
choking  occurs  and  a  shock  is  formed  within  the  combustor  and 
propagates  upstream. 

All  the  simulations  presented  are  carried  out  using  the  same 
sequence  of  events.  Initially  (t  e  [0  :  0.5  ms])  a  supersonic  flow  is 
established  in  the  engine.  Afterward,  the  fuel  is  injected  in  the 
combustor  with  a  constant  flow  rate  for  a  short  time  interval 
(f  e  [0  :  1.5  ms]).  Afterward  (t  >  1.5  ms)  no  fuel  is  injected.  As  an 
example,  two  different  values  of  the  fuel  mass  flow  rate  are 
considered  in  the  calculations  reported  in  Fig.  3.  At  the  low  fuel 
flow  rate  the  flow  remains  supersonic  throughout  the  scramjet;  in 
the  high  fuel  flow  rate  case,  thermal  choking  occurs  and  a  shock 
wave  moves  upstream  into  the  isolator  and  eventually  the  engine 
is  fully  unstarted. 

Interestingly  the  present  model  is  also  able  to  reproduce 
intermediate  conditions  as  illustrated  In  Fig.  4.  For  a  range  of  fuel 
flow  rates,  thermal  choking  occurs  (detected  by  the  presence  of 
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Fig.  3.  Time/space  plot  of  the  Mach  number  within  the  simplified  scramjet 
geometry.  Fuel  is  injected  with  a  constant  flow  rate  between  tc  [0.5  ms:  1.5  msj. 
Low  (top)  and  high  (bottom)  fuel  flow  rate  are  considered,  leading  to  a  stable 
operation  of  the  engine  and  to  unstart,  respectively.  The  dashed  black  line  is  the 
sonic  line  which  indicates  if  thermal  choking  occurred. 


subsonic  flow)  and  a  shock  is  formed  at  the  entrance  of  the 
combustor;  nevertheless  stable  burning  conditions  are  achieved;  a 
slight  increase  in  the  fuel  injection  results  in  the  shock  moving 
with  a  slow  velocity  upstream.  Once  the  fuel  injection  is  termi¬ 
nated  (t=1.5ms),  supersonic  conditions  are  again  established 
(see  Fig.  4)  in  both  cases  without  unstart. 

3.3.  Quantity  of  interest 

The  fundamental  feature  of  the  unstart  process  is  the  forma¬ 
tion  of  a  normal  shock  wave  traveling  upstream  in  the  isolator 
and  eventually  reaching  the  engine  inlet.  It  is  natural  to  consider 
the  location  and  motion  of  the  shock  as  the  primary  quantity  of 
interest.  Discontinuities  create  inherent  difficulties  for  numerical 
algorithms  and  in  the  present  approach  we  introduce  artificial 
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Fig.  4.  Time/space  plot  of  the  Mach  number  within  the  simplified  scramjet 
geometry.  Fuel  is  injected  with  a  constant  flow  rate  between  te[0.5  ms:  1.5  ms]. 
With  the  given  fuel  flow  rate  a  shock  initially  moves  upstream  in  the  isolator,  but 
supersonic  conditions  are  recovered  once  the  injection  is  terminated.  The  dashed 
black  line  is  the  sonic  line  which  indicates  if  thermal  choking  occurred. 


dissipation  in  the  form  of  a  low-order  discretization  in  the  vicinity 
of  shocks.  The  addition  of  dissipation  smears  the  discontinuity 
over  several  grid  cells  creating  the  additional  challenge  of  pre¬ 
cisely  defining  its  position.  In  this  work,  we  follow  [15]  by 
introducing  a  specific  algorithm  to  reconstruct  the  shock  location 
independently  of  the  computational  grid  used:  the  density  field  in 
the  vicinity  of  the  shock  is  reconstructed  using  a  hyperbolic 
tangent  fit.  At  any  given  time,  we  assume  a  representation  of 
the  density  profile  in  the  form: 


p(x)  = 


Pl  +  Pr  , 
2 


l^l-tanh 


(5) 


where  the  parameters  p^,  p^,  Xs  and  w  are  obtained  using  a  least 
square  minimization  procedure  starting  from  the  computed 
density. 


Table  1 

Convergence  rate  for  the  shock  tube  problem.  The  L2  norm  of  the  error  in  density  is 
defined  in  x  e  [3: 4]  and  t=2  and  p  is  defined  using  three  successively  refined  grids 
116],  Xs  is  the  shock  location  computed  using  the  hyperbolic  tangent  fitting,  to  be 
compared  to  the  exact  location  Xs=3. 5043116. 


Cells 

Ax 

At 

L2  X  10^ 

p 

100 

0.1000 

0.00200 

0.68020 

_ 

3.52890 

200 

0.0500 

0.00100 

0.34652 

- 

3.51720 

400 

0.0250 

0.00050 

0.18042 

1.00640 

3.51090 

800 

0.0125 

0.00025 

0.09992 

1.04499 

3.50760 

3.4.  Verification 

The  computational  algorithm  used  to  solve  Eq.  (1)  has  been 
applied  to  the  classic  time-dependent  shock  tube  problem.  Two 
chambers  (left,  L  and  right  R)  are  filled  with  stagnating  gas  and 
separated  by  a  diaphragm,  that  at  time  f=0  is  fractured.  The 
density  ratio  considered  is  Pi^/p^  =  8  and  the  pressure  ratio 
Pl/Pr  =  10. 

A  quantitative  comparison  between  the  exact  solution  and  the 
solutions  obtained  using  two  different  grids  (with  100  and  500 
cells)  is  given  in  Table  1  where  the  L2  norm  of  the  error  in  density 
in  the  region  around  the  shock  wave  (x£[3:4]  and  t=2)  is 
reported.  First  order  convergence  in  the  vicinity  of  the  disconti¬ 
nuity  is  expected  given  the  numerical  discretization  employed. 

The  shock  position  computed  using  the  hyperbolic  fit  is  also 
reported  in  Table  1,  showing  that  the  increased  grid  resolution 
corresponds  to  better  prediction  of  the  shock  location. 


3.5.  Calibration 

The  reduced  model  we  introduced  earlier  is  an  extreme 
simplification  of  the  complex  injection,  mixing  and  combustion 
processes  occurring  in  a  scramjet;  the  precise  distribution  of  heat 
release  has  a  strong  Impact  on  the  characteristics  of  unstart  and 
therefore,  it  is  important  to  calibrate  the  model  using  relevant 
experimental  data.  The  experimental  study  of  [17]  is  used  as  a 
reference.  The  geometry  is  a  simple,  constant  cross-section  duct; 
the  measurements  were  carried  out  for  different  equivalence 
ratios  and  inlet  Mach  numbers.  The  present  model  Is  calibrated 
using  one  set  of  data  (corresponding  to  Ma,= 2.5  and  (p  =  0.4)  and 
then  validated  using  two  different  operating  conditions  in  the 
following  section. 

The  objective  of  the  present  model  is  to  reproduce  the  relation 
between  the  heat  release  in  the  chamber  and  the  occurrence  of  a 
normal  shock  propagating  upstream  in  the  engine  isolator.  When 
comparing  to  measured  data,  the  initial  computed  location  of  the 
shock  always  corresponds  to  experimental  observations  (the 
initial  distance  from  the  actual  Injector);  this  allows  us  to  assess 
the  ability  of  the  model  to  reproduce  the  shock  propagation  speed 
which  ultimately  determines  whether  unstart  will  occur. 

The  shock  locations  measured  in  [17]  are  used  to  determine 
the  two  constants  (Q  and  C2)  present  in  the  heat  release  model. 
In  Fig.  5  the  results  obtained  using  the  calibrated  model  are 
compared  to  the  reference  data  used;  note  that  the  experimental 
uncertainty  is  mainly  associated  with  the  determination  of  the 
shock  location  extracted  from  Schlieren  images.  The  values  of  the 
constants  in  Eq.  (4  are  C]=4.21  and  C2=2.15;  these  correspond  to 
the  best  fit  of  the  experimental  shock  location.  Note  that  the 
values  used  In  [8]  were  C]  =  3  and  C2  =  l.  Computations  are  also 
repeated  assuming  +  5%  variation  in  the  two  constants  to  identify 
the  variability  associated  with  the  present  imprecise  calibration; 
the  results  are  reported  as  intervals  in  Fig.  5. 
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Fig.  5.  Calibration  results:  computations  of  the  shock  speed  for  inlet  Mach 
number  of  2.5  and  0  =  0.4  (see  [17]).  The  experimental  data  are  used  for 
calibrating  the  constants  in  the  combustion  model  (Eq.  (4)).  The  uncertainty  bars 
in  the  simulations  correspond  to  +  5%  change  in  the  model  calibration  constants. 
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Fig.  6.  Validation  results:  computations  of  the  shock  speed  for  inlet  Mach  number 
of  3.8  and  0  =  1.2  and  0  =  1.7  (see  [17])  using  the  previously  calibrated  model. 
The  uncertainty  bars  in  the  simulations  correspond  to  +  5%  change  in  the  model 
calibration  constants. 


3.6.  Validation 

The  calibrated  model  is  used  here  to  assess  its  capability  to 
reproduce  shock  propagation  speeds  at  various  equivalence 
ratios,  in  particular,  we  again  consider  the  experimental  data  set 
of  [17]  but  we  focus  on  higher  inflow  Mach  numbers  than  those 
corresponding  to  the  calibration  test.  As  before,  the  accreditation 
is  limited  to  the  evaluation  of  the  computed  shock  speed.  The 
results  for  two  values  of  the  equivalence  ratio  are  reported  in 
Fig.  6  for  an  inflow  Mach  number  of  3.5;  the  computations  are 
also  repeated  changing  the  constant  in  the  heat  release  model 
by  5%. 

The  computational  results  show  that  the  present  one-dimen¬ 
sional  model  represents  the  thermal  choking  and  the  subsequent 
normal  shock  formation  and  motion  in  a  reasonable  fashion,  at 
least  for  moderately  rich  mixtures  (</>  <  2.0).  Results  for  higher 
fuel  inflow  rates  could  not  be  obtained  because  of  numerical 
difficulties:  the  extreme  heat  release  leads  to  choking  immedi¬ 
ately  after  the  start  of  the  fuel  injection  phase.  A  very  strong 
shock  is  formed  and  this  leads  to  spurious  density  oscillations. 
We  conjecture  that  the  simple  heat  release  distribution  assumed 
in  the  chamber  (Eq.  (4))  is  not  credible  in  the  case  of  high  fuel 


injection  rate,  because  the  effectiveness  of  the  air/fuel  mixing  will 
likely  determine  an  upper  bound  on  the  effective  heat  release. 

4.  QMU:  quantification  of  margins  and  uncertainties 

The  first  step  in  developing  the  QMU  framework  is  to  define 
one  or  more  figures  of  merit  (metric)  to  characterize  successful 
operation  of  the  scramjet  engine.  It  is  important  that  the  chosen 
metrics  be  sensitive  to  all  the  important  aspects  of  the  prediction 
tools  (and  the  experiments)  used  in  the  evaluation.  In  the  present 
scenario  the  engine-generated  thrust  is  the  key  characteristic  of 
the  propulsion  system,  and  it  is  an  obvious  choice.  Additional 
quantities  could  also  be  selected,  such  as  the  structural  loads  in 
the  combustor,  but  for  sake  of  simplicity  only  one  metric  is 
considered  here. 

In  Fig.  7  a  hypothetical  thrust  signature  is  reported;  as  the  fuel 
flow  rate  is  increased  (or  equivalently  according  to  our  model  as 
the  equivalence  ratio  tp  is  increased),  combustion  takes  place  in 
the  engine.  When  a  sufficient  amount  of  heat  is  released  in  the 
mean  supersonic  stream,  the  pressure  increases  on  the  nozzle 
walls  and  leads  to  net  thrust.  As  the  fuel  rate  is  further  increased 
the  amount  of  thrust  can  be  adjusted.  In  Fig.  7  a  lower  bound 
corresponding  to  no  net  thrust  is  reported  together  with  an  upper 
bound  associated  with  thermal  choking  that  eventually  leads  to 
engine  unstart.  These  two  bounds  characterize  the  performance 
gate,  or  the  operating  region  in  Fig.  1 . 

The  confidence  ratio,  CR,  must  now  be  evaluated  with  respect 
to  this  critical  gate,  and  indicates  whether  the  system  is  operating 
safely  within  the  performance  gate;  therefore  we  need  to  char¬ 
acterize  both  the  margin  and  the  uncertainties. 

4.1.  Margins 

The  margin  is  loosely  defined  as  the  desired  distance  between 
the  performance  threshold  and  the  designed  operation  of  the 
vehicle.  As  mentioned  earlier,  from  a  design  perspective  the 
highest  thrust  is  achieved  when  operating  as  close  as  possible 
to  the  thermal  choking  limit. 

Detailed  experimental  investigations  have  been  carried  out  to 
identify  the  conditions  that  determine  choking  in  the  combustor; 
a  notable  example  is  the  work  of  Owens  et  al.  [18]  who  construct 
an  empirical  expression  for  the  critical  pressure  in  the  combustor 
in  the  form: 

Pcr=f(Mai,rhf,Ai/Ae,y)  (6) 


Fig.  7.  Performance  gate  for  the  scramjet  system. 
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where  Ma,  is  the  entrance  Mach  number,  ifif  is  the  fuel  flow  rate, 
A,  and  Ag  are  the  entrance  and  exit  cross-sectional  area  of  the 
combustor  (assumed  to  be  simply  a  diverging  duct)  and  y  is  the 
ratio  of  the  specific  heats. 

A  more  recent  study  [8]  approached  the  problem  from  a 
theoretical  perspective,  using  only  thermodynamic  laws  and 
quasi-lD  assumptions.  As  a  final  result  they  expressed  the  critical 
thrust  (corresponding  to  thermal  choking)  in  a  form  similar  to 
Eq.  (6). 

As  mentioned  earlier,  our  goal  Is  to  illustrate  the  development 
and  application  of  a  QMU  framework  and,  given  the  extreme 
simplification  of  the  computational  model.  It  Is  reasonable  to 
consider  a  rather  elementary  measure  of  the  thermal  choking 
limit.  We  follow  the  experimental  analysis  in  [19]  and  consider  an 
upper  bound  on  the  ratio  of  the  pressure  in  the  combustor  over 
the  inlet  pressure  as  the  indicator  for  unstart;  in  this  framework, 
Le  [19]  proposes  a  quadratic  dependency  of  the  bound  with 
respect  to  the  equivalence  ratio;  we  will  consider  the  same  model. 

The  complex  nature  of  the  supersonic  mixing  and  the  auto¬ 
ignition  process  is  not  well  understood  and  this  leads  to  uncer¬ 
tainties  In  the  determination  of  the  upper  bound  on  the  amount 
of  fuel  that  can  be  safely  injected  before  unstart  occurs.  In 
particular,  an  increase  in  fuel  flow  rate  can  lead  to  an  initial  large 
heat  release  that  triggers  the  inception  of  a  normal  shock  and  the 
formation  of  a  subsonic  region  in  the  combustor.  The  initial 
upstream  motion  of  this  shock  can  lead  to  boundary  layer 
separation;  the  resulting  highly  turbulent  flow  will  be  dominated 
by  large  vortical  structures  that  can  potentially  lead  to  a  reduction 
of  mixing  and  local  quenching,  e.g.  combustion  will  occur  only  in 
pockets  within  the  combustor.  The  resulting  reduction  in  the 
overall  heat  release  could  effectively  stop  the  upward  motion  of 
the  shock  and  lead  to  recovery  and  establishment  of  supersonic 
combustion  conditions.  In  other  words,  the  unsteady  response  of 
the  engine  to  an  increase  of  the  fuel  flow  rate  Is  critically 
controlled  by  the  presence  of  large  scale  structures  and  the 
system  might  continue  to  operate  normally  even  If  temporarily 
above  the  static  limit  of  heat  release.  The  upper  margin  of  Fig.  7  is 
therefore  itself  uncertain  as  it  is  controlled  by  time-dependent 
events  that  are  quite  stochastic  in  nature.  In  the  present  context, 
given  the  limited  information  available,  it  is  appropriate  to 
assume  that  the  gate  is  characterized  by  an  interval  [20],  as 
opposed  to  a  probability  distribution  function  [21].  Given  the 
limited  experimental  evidence  (for  example  [19,  Fig.  6])  we 
estimate  the  uncertainty  conservatively  as  +  8%  of  the  pressure 
ratio.  It  Is  worth  noting  that  the  experiments  in  [19]  only 
considered  a  combustor  geometry  characterized  by  a  diffuser 
angle  oicomb  =  2.5  ;  we  will  assume  that  the  diffuser  angle  In  the 
combustor  affects  this  bound  in  a  linear  fashion.  In  summary  the 
upper  bound  estimate  is  expressed  as 

^  =(a0^  +  b0  +  c)(a„m(i-l-5)  +  8%  (7) 

Pi 

where  a=18.6,  b  =  14.8  and  c=1.40. 

The  lower  bound  on  the  operating  conditions  is  related  to  the 
ability  to  generate  thrust,  overcoming  the  aerodynamic  resis¬ 
tance.  In  the  application  presented  we  will  Just  assume  a  low 
value  of  Pcr/Pi  as  being  the  no  net-thrust  limit. 

4.2.  Uncertainties 

Uncertainties  are  always  abundant  In  complex  engineering 
systems;  the  scramjet  is  no  exception.  Even  in  the  simplest 
schematic  representation  of  the  engine,  several  sources  of  uncer¬ 
tainty  must  be  accounted  for  as  they  can  dominate  the  overall 
performance.  In  the  present  context,  we  will  focus  on  aleatory 
uncertainties,  associated  with  imprecise  characterization  of  the 


Table  2 

Uncertain  quantities  considered  in  the  present  study  with  the  associated  mean 
and  assumed  variability.  Ma„  Tt,  and  y  characterize  the  inflow  conditions,  whereas 
C,  and  C2  the  uncertainty  in  the  heat  release  model  (Eq.  (4)). 


Quantity 

Type 

Mean 

Bounds {%) 

Ma,- 

Aleatory 

2.5 

±3 

Tt, 

Aleatory 

300 

±4 

7 

Aleatory 

1.4 

±5 

Cl 

Epistemic 

4.21 

±5 

C2 

Epistemic 

2.15 

±5 

various  parameters  affecting  the  engine  operation,  and  epistemic 
uncertainty  related  to  the  calibration  step  presented  earlier. 
In  particular,  we  consider  variability  in  flight  conditions  (inflow 
Mach  number  and  total  temperature),  the  thermodynamic  state 
of  the  incoming  air  (changes  in  }’  occurring  because  of  potential 
thermal  non-equilibrium  effects  at  the  bow  shock),  and  the 
injected  fuel  mass  flow  rate.  The  assumed  uncertainty  is  reported 
in  Table  2. 

The  present  computational  model,  on  the  other  hand,  is  a  very 
rudimentary  representation  of  the  rich  thermal/fluid  scenario 
within  the  engine.  We  are  neglecting  viscous  effects,  although  It 
is  well  known  that  the  interaction  of  shocks  with  boundary  layers 
plays  a  crucial  role  in  controlling  the  unstart  progression. 
In  addition,  the  entire  engine  is  assumed  to  operate  in  adiabatic 
conditions  and  with  fixed  geometry,  whereas  in  reality  the 
thermal  coupling  between  the  flow  and  the  structures  Is  very 
important.  We  plan  to  address  these  and  other  sources  of 
modeling  uncertainties  (epistemic  uncertainty)  in  the  future. 

4.3.  Confidence 

Evaluation  of  the  confidence  ratio  CR  requires  the  calculations 
of  measures  of  both  margins  and  uncertainties.  In  a  deterministic 
setting,  the  margin  is  simply  the  distance  between  the  perfor¬ 
mance  threshold  and  the  design  condition;  given  the  variability  in 
the  performance  estimate  and  the  uncertainty  in  the  unstart 
boundary,  we  must  distinguish  between  the  expected  margin 
related  to  the  best  estimate  of  the  bound  and  the  operating 
condition  and  the  effective  margin  which  accounts  for  the  pre¬ 
sence  of  uncertainty;  only  the  latter  is  of  practical  Importance. 
The  appropriate  definition  of  the  margin  depends  on  the  specific 
representation  of  the  uncertain  quantities. 

Sharp  et  al.  [20]  introduced  the  concept  of  full  and  interval 
QMU  to  characterize  the  confidence  in  the  two  extreme  cases 
when  either  probabilistic  descriptions  or  only  estimated  ranges  of 
the  quantities  of  interest  are  available.  Specifically,  we  define; 

•  Full  QMU:  The  analysis  is  based  on  complete  and  precise 
probabilistic  descriptions  of  all  the  quantity  of  interest.  In  this 
case,  the  effective  margin  relative  to  the  upper  boundary  can  be 
defined  as  (Fig.  8)  by  considering  that  for  both  the  design 
conditions  (subscript  D)  and  performance  threshold  (subscript 
T)  we  can  compute  the  mean  (p)  and  the  standard  deviation 
(cr)  of  the  corresponding  probability  distributions.  We  define; 

M  =  (pj-naT)-{po  +  n(^D):  U  =  n^a}  +  a1 

where  n  is  the  number  of  standard  deviation  values  that 
correspond  to  a  given  confidence  level.  Note  that  with  the 
above  definitions,  the  metrics  for  both  the  margin  M  and  the 
uncertainty  U  are  described  by  a  single  value.  Interestingly  in 
this  case  the  margin  can  be  negative  independent  of  the 
uncertainty  in  the  bounds  on  the  operability  limit  (inception 
of  unstart). 
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prediction 
(best  estimate) 


Fig.  8.  Upper  bound  of  the  performance  gate  for  the  scramjet  system  (cfr.  Fig.  7). 
The  expected  margin  is  defined  as  the  distance  between  the  performance  thresh¬ 
old  (subscript  T)  and  the  design  conditions  (subscript  D).  The  effective  margin 
accounts  for  the  uncertainty  in  the  evaluation  of  both  quantities.  Top:  Full 
QMU — probability  distribution  for  both  the  performance  threshold  and  the  design 
condition  estimate  are  available.  Bottom:  Interval  QMU — only  intervals  are 
available  for  both  quantities. 


•  Interval  QMU:  the  uncertainties  are  described  using  only 
ranges  (intervals)  without  referring  to  probabilistic  descrip¬ 
tions.  This  is  a  scenario  compatible  with  lack-of-information 
and  is  considered  in  the  analysis  that  follows.  In  this  case,  the 
uncertainty  in  the  lower  bound  of  the  performance  gate  is 
simply  given  by  Uj  while  the  uncertainty  in  the  upper  bound 
of  the  operating  region  is  given  by  Ud  (see  [20]  for  details  of 
how  to  handle  a  combination  of  intervals  and  probability 
distributions).  To  evaluate  the  margin  we  define  the  expected 


performance  threshold  (/tj-)  and  design  condition  (/tj,)  to  be 
the  mean  of  the  respective  intervals.  In  this  case  we  define: 

M  =  (llT  —  UT)—(ft[i  +  U[)):  U  =Ut  +  Uo 


5.  QMU  analysis  of  thermal  choking 

Simulations  are  carried  out  corresponding  to  the  nominal 
design  conditions  specified  in  Table  2.  One  of  the  problems  in 
propagating  the  identified  uncertainty  in  the  computational  code 
represented  by  Eq.  (1)  is  the  difficulty  in  representing  both 
probabilistic  (aleatory)  and  epistemic  uncertainty.  In  the  follow¬ 
ing  examples  we  considered  all  the  quantities  to  be  represented 
as  independent,  uniform  random  variables  and  we  used  Monte 
Carlo  sampling  to  characterize  the  output  variability  in  terms  of 
pressure  distribution  within  the  engine.  In  the  analysis  of  the 
results  we  will  consider  the  support  of  the  output  distributions  as 
an  interval  with  no  associated  probability  distribution.  It  is 
possible  to  retain  the  probabilistic  description  of  the  aleatory 
uncertainties  and  reflect  the  corresponding  distribution  in  the 
determination  of  the  margins;  one  possible  approach  is  to  use  a 
nested  sampling  strategy  [22].  It  is  important  to  point  out  that 
such  methodology  would  lead  to  an  evaluation  of  CR  that  would 
be  more  complex  that  what  afforded  by  either  interval  or 
full  QMU. 

Fig.  9  shows  the  pressure  time-signal  at  two  locations:  within 
the  nozzle  and  at  the  entrance  of  the  combustor.  The  solution  is 
highly  sensitive  to  the  input  variability,  especially  at  x=0. 

More  insight  into  the  output  uncertainty  is  presented  in 
Fig.  10,  where  the  coefficient  of  variation  (COV)  for  the  pressure 
is  reported  as  a  space/time  plot.  Note  that  in  the  present  interval 
analysis  the  COV  is  defined  as  the  midpoint  of  the  interval  over 
the  range.  Extremely  high  values  of  COV  (  r;  100%)  occur  close  to 
the  combustor  entrance.  These  are  representative  of  incipient 
unstart  conditions,  where  under  uncertainty  a  normal  shock 
moves  upstream. 


5.3.  Evaluating  the  confidence  ratio 

The  evaluation  of  the  confidence  ratio  CR  requires  the  calcula¬ 
tions  of  both  margins  and  uncertainties.  We  performed  Monte 
Carlo  simulations  to  account  for  the  uncertainties  introduced 
before,  and  we  used  a  fixed  geometiy  and  a  number  of  different 
fuel  injection  conditions.  The  objective  is  to  study  the  operating 
region  of  the  scramjet,  as  the  probability  of  engine  unstart 
increases  as  the  more  fuel  flow  rate  increases.  In  Fig.  1 1  the  results 
of  this  analysis  are  reported  in  terms  of  the  pressure  within  the 
combustor;  two  different  behaviors  emerged,  corresponding  to 
normal  operation  (supersonic  flow  throughout  the  system)  and  to 
the  occurrence  of  a  shock  in  the  combustion  chamber  and,  there¬ 
fore,  to  engine  unstart.  This  abrupt  change  of  behavior  is  comple¬ 
tely  determined  by  the  uncertainties  present,  since  the  nominal 
inputs  of  the  model  are  selected  according  to  safe  operating 
conditions.  From  Fig.  1 1  we  can  readily  evaluate  the  margin  and 
the  uncertainty  in  the  quantity  of  interest  for  each  value  of  the 
equivalence  ratio.  It  is  clear  that  the  uncertainty  in  the  unstart 
bound  also  grows  considerably  as  the  equivalence  ratio  is 
increased.  It  is  important  to  note  that  the  accuracy  of  the  predic¬ 
tions  during  an  unstart  event  is  questionable,  because  the  present 
model  ignores  viscous  interactions  that  will  dominate  the 
dynamics  of  the  upstream  moving  shock.  For  this  reason  in 
Fig.  1  l(bottom),  we  replaced  the  predicted  unstart  bound  (reported 
in  Fig.  n(top))  with  the  experimental  correlation  -  and  the 
reported  uncertainty  -  introduced  in  [18[.  This  results  is  a 
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Fig.  9.  Predictions  of  the  pressure  within  the  engine:  towards  the  nozzle  exit  (top) 
and  at  the  entrance  to  the  combustor  (bottom).  Bars  represent  95%  confidence 
bounds  obtained  performing  Monte  Carlo  simulations. 

somewhat  smaller  margin  at  low  equivalence  ratio,  but  a  con¬ 
siderably  smaller  uncertainty  at  higher  fuel  flow  rates. 

5.2.  QMU  for  system  design 

The  primary  objective  of  QMU  is  to  evaluate  the  confidence 
(and  reliability)  of  a  system  for  a  range  of  design  conditions.  In  the 
present  context,  we  are  not  targeting  a  particular  scramjet,  and  it 
is  therefore  more  useful  to  illustrate  how  a  QMU  framework  can 
be  used  to  evaluate  the  interplay  between  operability  limits  and 
design  choices.  The  thermal  choking  limit  is  strictly  a  thermo¬ 
dynamic  phenomenon:  the  deceleration  of  a  supersonic  stream  in 
the  presence  of  heat  addition  eventually  leads  to  the  formation  of 
a  normal  shock.  On  the  other  hand,  the  unstart  phenomena 
triggered  by  the  thermal  choking  is  determined  by  a  pure  fluid- 
dynamics  process:  the  upstream  motion  of  the  normal 
shock  within  the  isolator.  From  the  example  illustrated  before 
(see  for  example  Fig.  3(top))  it  is  clear  that  even  in  the  presence  of 


Fig.  10.  Space/time  plot  of  the  coefficient  of  variation  of  the  pressure  within  the 
scramjet  engine. 

thermal  choking  the  upstream  propagation  of  the  normal  shock 
can  be  avoided  by  considering  a  diffuser-shaped  combustor;  in 
fact,  in  theory  the  larger  the  angle  Oicomb  the  more  heat  release  can 
be  accommodated  in  the  combustor.  From  a  fluid  dynamics 
perspective  the  upstream  shock  propagation  is  inhibited  by  the 
adverse  pressure  gradient  environment. 

Our  objective  is  to  compute  the  QMU  confidence  ratio  for  a 
combination  of  equivalence  ratios  and  combustor  diffuser  angles, 
so  as  to  construct  the  safe  operation  region  for  a  generic  scramjet 
engine.  We  performed  computations  similar  to  those  illustrated 
earlier  where  the  fuel  is  injected  for  1ms  (0  #  0  e  [0.0004: 
0.0016]  ms)  at  a  constant  equivalence  ratio.  We  define  unstart  to 
occur  when  a  shock  is  detected  at  a  certain  location  in  the  isolator 
(x  =  -0.4  m)  within  a  time  window  of  4  ms  (t  e  [0: 0.004]  ms).  The 
computations  are  repeated  for  (p  e[0-A  .6]  and  iXcomb  e  [0: 12]. 

The  results  are  presented  in  terms  of  a  confidence  ratio  map, 
reported  in  Fig.  12.  The  points  where  Ck>  1  correspond  to  safe 
operation  conditions.  It  is  worth  noting  that  two  artificial  bound¬ 
aries  have  been  added  to  the  plot  in  Fig.  12.  The  right  boundary, 
namely  the  flow  separation  limit,  is  related  to  the  eventual  flow 
separation  induced  by  large  diffusion  angles  in  the  combustor. 
The  present  model  is  unable  to  capture  this  phenomenon,  but  it  is 
expected  that  after  a  certain  angle  the  interaction  of  the  shock  in 
the  combustor  and  the  boundary  layers  will  create  a  large 
separation  that  eventually  will  lead  to  loss  of  thrust;  we  assumed 
for  simplicity  that  this  bound  corresponds  to  a  given  diffuser 
angle.  The  second  boundary  (represented  by  the  lower  curve  in 
Fig.  12)  indicates  the  minimum  value  of  the  thrust  that  is  required 
for  acceptable  vehicle  performance.  This  bound  can  only  be 
evaluated  if  additional  information  on  the  system  configuration 
is  given;  here  we  assumed  that  the  thrust  is  directly  proportional 
to  the  injected  fuel  (</>)  and  that  the  diffuser  angle  in  the 
combustor  contributes  to  create  thrust.  In  the  present  analysis 
we  assumed  that  these  two  bounds  are  known  as  they  are  not 
directly  related  to  the  thermal  choking  phenomenon. 

The  CR  factor  is  evaluated  by  performing  Monte  Carlo  simula¬ 
tions  at  each  point  in  the  (p-o-comb  plane;  in  practical  terms 
only  8x8  analyses  (each  one  consisting  of  1000  independent 
realizations  aimed  at  characterizing  the  uncertainty  described  in 
Table  2)  were  carried  out  and  the  results  interpolated  to  construct 
the  curves  in  Fig.  12.  The  results  show  several  interesting  features. 
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Fig.  11.  Evaluation  of  unstart  margin  for  a  range  of  fuel  equivalence  ratios.  The 
unstart  boundary  is  evaluated  using  the  current  computational  model  (top)  and  an 
experimental  correlation  (bottom)  [18], 

the  two  most  important  observations  from  a  scramjet  operation 
perspective  are: 

1.  For  a  straight  combustor  geometry  {of.comb  =  0)  the  operating 
region  is  very  narrow  in  terms  of  the  injected  fuel  equivalence 
ratio.  Increasing  the  demand  on  the  reliability  (demanding  a 
larger  confidence  ratio)  leads  to  an  unfeasible  design. 

2.  The  safe  operating  region  expands  with  oCcomb-  physically  this  is 
due  to  the  increased  ability  to  stabilize  a  normal  shock  in  the 
combustor  resulting  from  the  diffuser  angle  in  the  combustor. 
For  large  diffuser  angles  an  increased  reliability  demand  does 
not  shrink  the  operability  range  considerably. 

Finally,  it  is  useful  to  note  that  the  presence  of  two  other 
boundaries  in  Fig.  12  naturally  leads  to  the  need  to  define  a  more 
comprehensive  model  of  the  scramjet  fiow/thermal  environment 
than  the  one  employed  in  this  study.  The  lower  bound  on  the 
vehicle  thrust  can  be  defined  according  to  typical  design  con¬ 
siderations  [12];  on  the  other  hand,  the  bound  on  the  right  -  the 


Fig.  12.  Unstart  boundaries  due  to  thermal  choking  for  two  required  confidence 
ratios  as  a  function  of  the  equivalence  ratio  and  the  diffuser  angle  in  the 
combustor.  The  grey  area  corresponds  to  the  safe  operation  region  (flj  in  Fig.  1). 
The  minimum  thrust  and  the  flow  separation  bounds  are  assumed. 

flow  separation  limit  -  requires  a  considerably  more  sophisti¬ 
cated  model  involving  the  solution  of  two-  or  three-dimensional 
reacting  flow  equations.  In  addition  the  QMU  analysis  also 
requires  the  determination  of  a  theoretical/experimental  limit 
for  the  corresponding  bound,  which  in  the  case  of  flow  separation 
can  be  extremely  challenging. 

From  a  computational  perspective,  the  construction  of  a  map 
such  as  that  shown  in  Fig.  12  is  a  formidable  undertaking.  64,000 
independent  unsteady  simulations  were  carried  out;  given  the 
extreme  simplicity  of  the  present  reduced-order  model  each 
computation  is  completed  in  less  than  10  s  and,  therefore,  all 
the  simulations  were  completed  in  about  a  day  on  an  8-cores 
desktop.  Each  solution  corresponds  to  a  point  in  the  4>-a^omb 
plane  and  provides  one  evaluation  of  the  (possible)  presence  of  a 
shock  traveling  upstream  towards  the  scramjet  inlet. 

In  Fig.  12  an  additional  line  corresponding  to  Ck=2  is  reported. 
This  line  indicates  the  operating  conditions  that  lead  to  an 
increased  performance  margin,  or  equivalently  that  are  less 
dominated  by  uncertainties.  In  the  present  context  the  CR 
represents  a  metric  that  is  equivalent  to  an  engineering  safety 
margin  and  can  be  used  as  a  free  parameter  to  explore  different 
configurations  and  design  choices.  An  additional  interpretation  of 
the  difference  between  the  CR=2  and  the  Ck=l  lines  is  the 
performance  gain  that  can  be  obtained  by  reducing  the  uncer¬ 
tainty  by  50%;  the  results  in  Fig.  12  show  that  for  low  combustor 
angles  the  increase  in  the  amount  of  fuel  allowed  without 
unstarting  the  engine  is  considerable,  whereas  only  a  limited 
change  is  observed  for  large  value  of  o-comb- 


6.  Conclusions 

The  QMU  analysis  presented  here  is  the  first  step  towards  a 
comprehensive  analysis  of  failure  modes  in  a  scramjet  engine  using 
computation.  We  define  a  simplified  simulation  framework  capable 
of  reproducing  qualitatively  thermal  choking  and  the  eventual 
engine  unstart.  This  model  was  verified  and  then  calibrated  against 
experimental  data  corresponding  to  a  simple  straight  duct  config¬ 
uration.  A  validation  step  was  also  carried  out  to  assess  the  accuracy 
of  the  simulations.  In  a  second  phase,  uncertainty  in  both  the  inlet 
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conditions  and  the  calibration  constants  in  the  heat  release  model 
were  evaluated  to  Identify  their  impact  on  the  unstart  scenario. 
Finally  the  performance  threshold  corresponding  to  unstart  trig¬ 
gered  by  thermal  choking  was  identified  and  analyzed. 

Future  work  will  focus  on  the  introduction  of  a  more  compre¬ 
hensive  computational  model  to  represent  the  reactive  flow  in  the 
scramjet.  The  inclusion  of  viscous  effects  is  clearly  of  primary 
importance  in  the  characterization  of  the  unstart  boundaries, 
together  with  a  better  representation  of  the  mixing,  combustion 
and  heat  release  In  the  combustion  chamber.  In  addition,  we 
intend  to  pursue  a  more  complete  characterization  of  the  mar¬ 
gins,  by  introducing  multiple  gates  that  represent  in  more  detail 
the  sequence  of  physical  processes  that  contribute  to  determine 
the  performance  thresholds.  The  extension  of  present  QMU 
framework  to  include  multiple  gates  Is  a  very  important  and 
challenging  endeavor. 
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